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We report an observation of the decay B° — > D*~t + v t in a data sample containing 535 x 10 6 
BB pairs collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. We 
find a signal with a significance of 5.2<r and measure the branching fraction B(B° — > D*~t + v t ) = 
(2.02t.o'1yr(stat)±0.37(syst))%. This is the first observation of an exclusive B decay with a b — > crv T 
transition. 

PACS numbers: 13.20.He, 14.40.Nd 



B meson decays with b — > ctv t transitions can pro- 
vide important constraints on the Standard Model (SM) 
and its extensions. Due to the large mass of the lep- 
ton in the final state these decays are sensitive probes of 
models with extended Higgs sectors [l| and provide ob- 
servables sensitive to new physics, such as polarizations, 
which cannot be accessed in other semileptonic decays. 

Multiple neutrinos in the final states make the search 
for semi-tauonic B decays very challenging and hence 
there is little experimental information about these pro- 
cesses. So far, results are limited to inclusive and 
semi-inclusive measurements by LEP experiments [H 
which measure an average branching fraction of B(b — > 
tv t X) = (2.48 ± 0.26)% SM calculations predict 
branching fractions for B — > D*t + v t around 1.4% with 
uncertainties arising mainly from assumptions about 
form- factors 

In this paper we present the first observation of 
B° — > D*~t + v t Q decay using a data sample contain- 
ing 535 x 10 6 BB pairs that were collected with the Belle 
detector at the KEKB asymmetric-energy e + e~ (3.5 on 
8 GeV) collider Q operating at the T(45) resonance 
{^fs = 10.58 GeV). The Belle detector is a large-solid- 
angle magnetic spectrometer consisting of a silicon ver- 
tex detector, a 50-layer central drift chamber, a system 
of aerogel Chercnkov counters, time-of-flight scintillation 
counters and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to identify K® mesons and muons. A detailed descrip- 
tion of the detector can be found in Rcf. T]. We use 
Monte Carlo (MC) simulations to estimate signal effi- 
ciencies and background contributions. Large samples 
of the signal B° — > D*~t + v t decays are generated with 



the EvtGcn package 0] using the ISGW2 model @. Ra- 
diative effects are modeled by the PHOTOS code [Toj ] . 
MC samples equivalent to about twice the accumulated 
data are used to evaluate the background from BB and 
continuum qq (q = u, d, s, c) events. 

B decays to multi-neutrino final states can be observed 
at B-factories via the recoil of the accompanying B me- 
son (Stag) [HI- Reconstruction of the Btag strongly sup- 
presses the combinatorial and continuum backgrounds 
and provides kinematical constraints on the signal me- 
son (B s i g ). In this study we take advantage of the clean 
signature provided by the D* meson occurring on the 
signal side and reconstruct the -Btag "inclusively" from 
all the particles that remain after selecting candidates 
for B s i g daughters. We apply the analysis to B s i g de- 
cay chains that combine a high reconstruction efficiency 
with a low background level. The D*~ mesons are re- 
constructed in the D*~ — > D°tt~ decay channel. The 
D°'s are reconstructed in the K + ir~ and K + ir~ir° final 
states. The r + — * e + i> e D T and t + — > Tt + D T modes are 
used to reconstruct r lepton candidates. We do not in- 
clude the t + — » ii + v^D T mode because in the relevant 
momentum range the muon identification is inefficient. 
The t + — > ir + i> T channel has higher combinatorial back- 
ground than the purely leptonic mode, but the single neu- 
trino in r decay provides better kinematical constraints. 
For this mode we analyze only the D° — > K + ir~ decay. 

We select charged tracks with impact parameters that 
are consistent with an origin at the beam spot, and hav- 
ing momenta above 50 MeV/c in the laboratory frame. 
Muons, electrons, charged pions, kaons and (anti)protons 
are identified using information from particle identifica- 
tion subsystems. The electrons from signal decays are se- 
lected with an efficiency greater than 90% and a misiden- 
tification rate below 0.2%. The momenta of particles 
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identified as electrons are corrected for brcmsstrahlung 
by adding photons within a 50 mrad cone along the 
trajectory. The ir° candidates are reconstructed from 
photon pairs having invariant mass in the range 118 
MeV/c 2 < M 77 <150 MeV/c 2 . From candidates that 
share a common 7, we select the tt° with the smallest x 2 
value from a mass-constrained fit. To reduce the combi- 
natorial background, we require photons from the tt° to 
have energies above 60 MeV -120 MeV, depending on the 
photon's polar angle. Photons that do not come from a 
7T° and exceed a polar-angle dependent energy threshold 
(100 MeV - 200 MeV) are included in the B tag recon- 
struction. 

We reconstruct the signal decay by selecting combina- 
tions of a D*~ meson and an electron or a pion candidate 
with opposite charge. We accept D° candidates with in- 
variant masses in a 5cr window around the nominal PDG 
Q value. D*~ candidates are accepted if the mass dif- 
ference Mq. — Mpo is in a 3cr window around the PDG 
value. In order to reduce background from incorrectly re- 
constructed tracks, we impose tighter impact parameter 
requirements on the e and tt candidates from t decay. 

Once a B s j g candidate is found, the remaining particles 
are used to reconstruct the B tag decay. The consistency 
of a -B tag candidate with a .B-mcson decay is checked us- 
ing the beam-energy constrained mass and the energy 

difference variables: M tag = ^ E^ cam - p 2 ag , p tag = 

YVpi, an d AE'tag = E tag — ^beamj E tag = Y^i^ii 

where Ebcam is the beam energy and Pi and Ei de- 
note the momentum vector and energy of the i'th par- 
ticle in the T(45) rest frame. The summation is over 
all particles that arc not assigned to -B s j g and satisfy 
the selection criteria described above. We require that 
events have at least one (D*~e + /ir + ) pair and that M tag 
and S tag satisfy M tag >5.2 GcV/c 2 and |A.E tag | <0.6 
GeV. To improve the quality of the -Btag reconstruction, 
we impose the following requirements: zero total event 
charge, no /.t^ and no additional in the event, zero 
net proton/antiproton number, residual energy in the 
ECL (i.e. the sum of energies of clusters that do not 
fulfill the requirements imposed on photons) less than 
0.35 GeV and number of neutral particles on the tag- 
ging side N^o + 7V 7 <5. These criteria, which we refer to 
as "the B t ag-selection" , reject events in which some par- 
ticles were undetected and suppress events with a large 
number of spurious showers. In order to validate the B tag 
simulation and reconstruction, we use a control sample 
of events, where the -B s i g decays to D*~n + (followed by 
D*~ — ► 5°7r~, D° — ► K + tt~) which allows us to se- 
lect a BB sample with a purity of 96% and with B slg 
and £>tag daughters properly assigned to the parent par- 
ticles. Figure 1 shows the M tag and AE tag distributions 
of the control sample for data and the MC simulation 
scaled to the integrated luminosity in data. The events 
satisfy the B tag -selection criteria and are in the —0.25 



GeV< A^tag < 0.05 GeV (for Fig. Ufa)) and M tag >5.27 
GeV/c 2 (for Fig. QJb)) windows. The good agreement 
of the shapes and of the absolute normalization demon- 
strates the validity of the MC-simulations for B tag de- 
cays. Based on this study we constrain all further anal- 
ysis to the region -0.25 GeV< AE tag <0.05 GeV. With 
this requirement about 80% of the events are contained 
in the range M tag > 5.27 GeV/c 2 . 




M taQ [GeV/c 2 ] AE [GeV] 



FIG. 1: Mtag and A-E tag distributions for B° — > D*~n + con- 
trol sample from data (points with error bars) and MC (his- 
tograms) . 

The procedure described above, when applied to events 
with (D*~e + ) pairs selects a relatively clean sample of 
semileptonic B decays with the dominant non-signal con- 
tribution from the B° — » D*~e + v e mode. Combinato- 
rial background from hadronic -B-decays dominates in 
the t + — > n + v T mode. The background suppression 
exploits observables that characterize the signal decay: 
missing energy E mis = E hcam - E D * - -B e/7r ; visible 
energy E v i a , i.e. the sum of the energies of all parti- 
cles in the event; the square of missing mass -M mis = 
-^mis — (Psig — PD' — Pc/tt) 2 an d the effective mass of the 
(tv t ) pair, = (E hcam - E D . f - (p sig - p^,, ) 2 where 
Psig = ~ Ptag- The most powerful variable for separating 
signal and background is obtained by combining E^ s 
and (D*e/w) pair momentum: X m i s = (E m { s — \pd* + 

Pe/irD/y^bcam m B° wnere m B° is the B° mass. The 
X mls variable is closely related to the missing mass in the 
B sig decay and does not depend on B tag reconstruction. 
It lies in the range [— 1 , 1] for events with zero missing 
mass (e.g. with a single neutrino) and takes larger val- 
ues if there are multiple neutrinos. The MC distributions 
of .Xmis and -E v i s for signal and background events after 
Btag-selection for the t — ► evv mode are shown in Fig. [2] 
The relative normalizations of the main background cat- 
egories, B° — ► D*~e + i/ e , B — > D**e + v e . other B decays 
and qq continuum, are determined from the data using 
looser selection criteria and verified using the sideband 
regions of the data sample that passed the final signal 
selection. 

We optimize selection criteria using MC samples for 
signal and backgrounds, separately for decay chains with 
t — > evv and with r — > nv. In the first case we 
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FIG. 2: Xmis and i? v is distributions (normalized to unity) 
after the Stag-selection for signal (blank) and background 
(shaded) for the r — * evv mode in the region M tag > 5.27 
GeV/c 2 . The background components, from top to bottom: 
B — ► D*~e + v e , B — > D**e + v e , and other B decays. The 
contribution from qrg-continuum is negligible. 



require A mis >2.75, 1.9 GeV< E mis <2.6 GeV and 
E viB <8.3 GeV. We also reject events with a small dif- 
ference between and M^ is to suppress background 
from hadronic B decays where a genuine D* meson 
is combined with a soft secondary e . Decays in the 
t — > 7tj/ mode are selected by requiring X m - ls >1.5, 
M^f — A/^is — m 2 T + m% >0 (m T and m, denote the 
masses of the r and charged 7T, respectively), E vis <8.3 
GeV, the energy of the tt + from the (D*~n + ) pair greater 
than 0.6 GeV, no K° L in the event and less than four 
tracks that do not satisfy the requirements imposed on 
the impact parameters. The second requirement is equiv- 
alent to the condition | cos0„ lt/2 | <1, where 6 VlV2 denotes 
the angle between the two neutrinos in the (r + v T ) rest 
frame. The last three criteria reduce combinatorial back- 
ground from low momentum pions and background from 
hadronic B — > D*~ K° L + X and B — > D*~nn+X decays. 
The above requirements result in flat M tag distributions 
for most background components, while the signal dis- 
tribution remains unchanged. This allows us to use the 
M tag variable to extract the signal. 

The M tag distribution of the signal is described using 
a Crystal Ball (CB) lineshape function [l^]. The shape 
parameters of the CB-function are determined from un- 
binned maximum likelihood fits to the combined MC sig- 
nal samples. All the fits are performed in the range 
Mtag > 5.2 GeV/c 2 . The backgrounds are modeled as 
the sum of a combinatorial component using a parame- 
terization introduced by ARGUS (ARGUS-function) [l]| 
and a peaking background described by the CB-function 
with shape parameters fixed from the fit to the signal 
MC. The main source of the peaking background is the 
semileptonic decay B° — ► D*~e + v e . Cross-feed events 
from signal decays followed by r decays to other modes 
are negligible in the t — > evv mode, but give significant 
contributions to the r — > nv mode. About half of the 
cross-feed comes from r — > pv decay. We parameterize 
the M tag distribution of cross-feed events as a sum of 



CB and ARGUS functions with shape parameters fixed 
from fits to the signal and combinatorial background as 
described above. The component described by the CB- 
function is treated as a part of the signal. The efficiencies 
of signal reconstruction and the expected combinatorial 
and peaking backgrounds are given in Table U 

The selection criteria established in the MC studies are 
applied to the data. The resulting Mt ag distribution for 
data in all three decay chains is shown in Fig. 03 The 
overlaid histogram represents the expected background, 
scaled to the data luminosity. A clear excess over back- 
ground can be observed. 




M tag [GeV/c?] 

FIG. 3: Mt a g distribution for the combined data sample. 
The histogram represents expected background scaled to the 
data luminosity. The solid curve shows the result of the fit. 
The dotted and dashed curves indicate respectively the fitted 
background and the component described by the ARGUS- 
function. 

We extract signal yields by fitting the M tag distribu- 
tions to the sum of the expected signal and background 
distributions using the following likelihood function: 

N 

C = e -(N s +N p+ N b ) -Q [(7Vs + Np)Ps{Xi) + Nb P b ( Xi )}, (1) 
i=l 

where Xi is the M tag in the i'th event and TV is the to- 
tal number of events in the data. P s (P&) denotes the 
signal (background) probability density function (PDF), 
which is parameterized as a CB (ARGUS )-function with 
shape parameters determined from fits to MC samples 
and N s , N b , and N p are the numbers of signal, com- 
binatorial background and peaking background respec- 
tively. N s and N b are free parameters of the fit, while 
N p is fixed to the value obtained from fits to MC sam- 
ples and scaled to the data luminosity (N p is set to 
zero for the r — ► iris mode). The fits are performed 
both for the three decay chains separately and for all 
chains combined with a constraint to a common value of 
B(B° -»■ D*-t+v t ). The fit results arc included in Ta- 
ble |U The total number of signal events is 60^^ with a 
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statistical significance of 6.7a. The significance is defined 
as £ = \f— 21n(£o / Anax) , where £ max and Cq denote 
the maximum likelihood value and the likelihood value 
for the zero signal hypothesis. The fitted signal yield 
is used to calculate the branching fraction for the decay 
B° — > D*~t + v t using the following formula, which as- 
sumes equal fractions of charged and neutral B mesons 
produced in T(45) decays: B = N S /(N BB x J2ij e ijBij)i 
where N BB is the number of BB pairs, e»j denotes 
the reconstruction efficiency of the specific decay chain 
and Bij is the product of intermediate branching frac- 
tions B(D*- -» D°tt-) x B(D° -»i)x B(t+ -> j). 
All the intermediate branching fractions are set to the 
PDG values Q. The branching fraction obtained is 
B(B° £*-t+z/ t ) = (2.02±° f 7 (stat))%. 

As a consistency check we also examine the distribu- 
tions used in the signal selection, applying all require- 
ments except those that are related to the considered 
variable. In all cases the distributions are well repro- 
duced by the sum of signal and background components 
with normalizations fixed from the fits to the M tag distri- 
bution. We also use the M^ is and cos9 Vl „ 2 (for r — > ttv 
mode) variables to extract the signal yield. We per- 
form fits to distributions of these variables in the region 
Mtag > 5.27 GeV/c 2 and obtain branching fractions in 
the range 1.83% - 2.05% and in agreement with the re- 
sults from the M tag fit. 

We consider the following sources of systematic un- 
certainties in the branching fraction determination. The 
systematic error on N BB is 1.3%. The systematic un- 
certainties in the signal yield arise from uncertainties 
in the signal and background shape and peaking back- 
ground. The systematic error due to the statistical un- 
certainties in the CB shape is 2.8%. The CB parameters 
obtained from MC-samples are, within statistical errors, 
consistent with those extracted from fits to the control 
sample in data. Therefore we do not introduce addi- 
tional uncertainties due to imperfect signal shape mod- 
eling. The systematic errors due to the parameteriza- 
tion of the combinatorial background are evaluated by 
changing the ARGUS-shape parameters by ±1ct. Fits 
with the shape parameters allowed to float provide con- 
sistent results within statistical uncertainties. The total 
systematic uncertainty due to the combinatorial back- 
ground parameterization is i^ 7 1 %. The systematic error 
due to the peaking background is evaluated for each chan- 
nel and amounts to ^44% for combined modes, which 
is dominated by MC statistics. The effective efficiency 
CijBij includes uncertainties in determination of the 
efficiencies for _B ta g reconstruction, (D*~ e + /ir + ) pair se- 
lection and signal selection. The uncertainty in Stag re- 
construction is taken as the statistical error in the Stag 
efficiency evaluated from the data control sample (tagged 
with B° — > D*~ir + decay) and is 10.9%. The systematic 
error on the determination of (D*~ e + /t: + ) pair selec- 



tion efficiency comes from systematic uncertainties in the 
tracking efficiency, neutral reconstruction efficiency and 
particle identification and is in the range 7.9%-10.7% de- 
pending on the decay chain. Systematic uncertainties 
in the signal selection efficiency are determined by com- 
paring MC and data distributions in the variables used 
for signal selection. The uncertainties due to the partial 
branching ratios are taken from the errors quoted in the 
PDG . All of the above sources of systematic uncer- 
tainties are combined together taking into account cor- 
relations between different decay chains. The combined 
systematic uncertainty is 18.5%. 

We include the effect of systematic uncertainty in the 
signal yield on the significance of the observed signal by 
convolving the likelihood function from the fit with a 
Gaussian systematic error distribution. The significance 
of the observed signal after including systematic uncer- 
tainties is 5.2<7. 

In conclusion, in a sample of 535xl0 6 BB pairs we 
observe a signal of 6OI44 events for the decay B° — > 
D*~t + v t with a significance of 5.2. This is the first ob- 
servation of an exclusive B decay with the b — > ctv t 
transition. The measured branching fraction: B(B° — > 
D*-t+v t ) = (2.02±gJ? (stat) ±0.37 (syst))% is consistent 
within experimental uncertainties with SM expectations 
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